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Abstract. Transverse energy, charged particle pseudorapidity distributions and photon transverse momen- 
tum spectra have been studied as a function of the number of participants {Npart) and the number of binary 
nucleon-nucleon collisions {Ncoii) in 158-yl GeV Pb+Pb collisions over a wide impact parameter range. A 
scaling of the transverse energy pseudorapidity density at midrapidity as ~ N^'^^^^'^'^^ and ~ -^coff ^° 

observed. For the charged particle pseudorapidity density at midrapidity we find a scaling as ~ N^'^Jf^^''^'^ 

and ~ N^cM^^^^^ ■ This faster than linear scaling with Npart indicates a violation of the naive Wounded 
Nucleon Model. 
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1 Introduction 

Heavy-ion collisions at ultrarelativistic energies probe nu- 
clear matter at high temperatures and densities. A major 
goal of these studies is the search for a deconfined phase 
of nuclear matter. A necessary condition to reach such a 
phase transition is local equilibration as might be achiev- 
able through rescattering of the produced particles. Since 
the amount of rescattering should increase with the size 
of the reaction system, it is of interest to study these re- 
actions over a wide range of centralities. 

For hard processes, where cross sections are small, the 
naive expectation is a scaling of the particle yields with the 
number of binary collisions. Experimentally, the scaling of 
cross sections with target mass in p-l-A collisions was stud- 
ied and the scaliiig was observed to be even stronger than 
this expectation M. This was later attributed to multiple 
parton scattering in the initial state From the same 
experiment it was also seen that particle production at in- 
termediate pt shows a much weaker increase with target 
mass. 

The gross features of particle production in nucleon- 
nucleus collisions and reactions of light nuclei are well de- 
scribed in the framework of the Wounded Nucleon Model 
Q. In this model the transverse energy and particle pro- 
duction in p-|-A and A-l-A reactions is calculated by as- 
suming a constant contribution from each participating 
nucleon. This kind of scaling has also been observed by the 
WA80 collaboration in reactions of ^^O and '^^S projectiles 
with various targets where dEx/drilmax was found to de- 
pend approximately linearly on the average total number 
of participants ||] . 

While a scaling with the number of collisions arises 
naturally in a picture of a superposition of nucleon-nucleon 
collisions, with a possible modification by initial state ef- 
fects, the Wounded Nucleon Model or participant scal- 
ing is more naturally related to a system with strong fi- 
nal state rescattering, where the incoming particles lose 
their memory and every participant contributes a similar 
amount of energy to particle production. The scaling be- 
havior of particle production may therefore carry impor- 
tant information on the reaction dynamics. Various exper- 
imental signatures in heavy ion reactions require a com- 
parison of observables for different system sizes. Therefore 
it is important to have a good understanding of these ba- 
sic scaling properties. The scaling behavior can also be 
used as a valuable test for models of particle production 
in heavy ion reactions (see e.g. |^). 

Furthermore, several observables in heavy ion reac- 
tions seem to show qualitative changes once a certain sys- 



tem size is reached. Strangeness production is enhanced 
in S-l-S reactions compared to P+P, but seems to saturate 
for even larger nuclei (see e.g. 0). Recent results from 
the WA98 experiment B show a significant change of the 
shape of the tt*^ pt spectrum in peripheral Pb-t-Pb colli- 
sions compared to p+p data. The shape, however, remains 
unchanged in the range of semi-central Pb-f Pb collisions 
with about 50 participating nucleons up to central reac- 
tions. 

The NA50 collaboration has observed an anomalously 
suppressed J/ip yield in central Pb+Pb collisions in con- 
trast to peripheral reactions |]l0|| where the suppression of 
the J/ip yield can be explained by absorption in nuclear 
matter. This anomalous J/ijj suppression provides an ad- 
ditional incentive to study the scaling behavior of particle 
production with the number of participants. Most models 
based on J/ip absorption by hadronic comovers assume 
the comover density to scale linearly with the number of 
participants and are then not able to fit the anomalous 
suppression pT| , p^ . Only if the hadronic comover density 
scales substantially faster than linearly with the number 
of participants it is possible to obtain reasonable fits to 
the anomalous suppression. It is therefore of interest to 
study in detail the centrality dependence of particle pro- 
duction and investigate its scaling properties with respect 
to the number of participants or collisions. 



2 Experiment and Data Analysis 

The CERN experiment WA98 is a general-purpose ap- 
paratus which consists of large acceptance photon and 
hadron spectrometers together with several other large 
acceptance devices which allow to measure various global 
variables on an event-by-event basis. The experiment took 
data with the 158-^ GeV ^ospb beams from the CERN 
SPS in 1994, 1995, and 1996. The layout of the WA98 
experiment as it existed during the final WA98 run pe- 
riod in 1996 is shown in Fig. [|. The data presented here 
were taken during the 1996 lead beamtime. The trans- 
verse energy and charged particle distributions shown in 
this paper were measured with the magnetic field of the 
Goliath magnet turned off. The minimum bias cross sec- 
tion for this configuration was amb = (6260 ± 280) mb. 
The error of the minimum bias cross section relates to the 
uncertainty of the target thickness and the uncertainty 
in the subtraction of the contribution from interactions 
outside the target. The contribution of these interactions 
was determined in special target-out runs. A 20% uncer- 
tainty of this contribution was assumed in the calculation 
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Fig. 1. The WA98 experimental setup. 



of the uncertainty of the minimum bias cross section. The 
WA98 experiment took most of the data with the Gohath 
magnet turned on. Since the analysis of the photon and 
neutral pion scaling requires high statistics the respective 
spectra used here were measured in a field-on configura- 
tion. The minimum bias cross section for this data set was 
cr,„fc = (6440 ± 300) mb. 

The Zero Degree Calorimeter is located 30 m down- 
stream of the target and measures the total energy of all 
particles within an angle < 0.3° relative to the beam 
axis in the laboratory system. The MIRAC calorimeter is 
placed 24 m downstream of the target It consists of 
a hadronic and an electromagnetic section and covers the 
pseudorapidity interval 3.5 < 77 < 5.5. MIRAC plays the 
central role in the WA98 minimum bias trigger where the 
measured Et is required to be above a minimum thresh- 
old. The systematic errors of dET/drj at midrapidity are 
dominated by the correction for the differences in the re- 
sponse of the MIRAC to hadronic and electromagnetic 
showers and to the extrapolation of the distribution of 
dET/drj to midrapidity. These combine to give an overall 
systematic uncertainty of « 20% in the absolute result for 
dEx / drjlmax ■ The centrality dependent part of this uncer- 
tainty is much smaller and estimated to be approximately 
5% only. The correction of dEx/d-q due to interactions 



outside the target for peripheral Pb-|-Pb reactions is less 
than 2%. 

The charged particle multiplicity is measured with a 
circular Silicon Pad Multiplicity Detector (SPMD) located 
32.8 cm downstream of the target [Q. It consists of four 
quadrants each produced from a 300 /im thick silicon wafer. 
This detector provides full azimuthal coverage of the pseu- 
dorapidity region 2.35 < 77 < 3.75 with 180 ^-bins and 22 
77-bins. The pad size increases radially to provide an ap- 
proximately uniform pseudorapidity coverage. In central 
Pb+Pb collisions the probability that a pad is hit by two 
or more particles is not negligible. Therefore the multiplic- 
ity in an 77-ring is determined from the sum of the mea- 
sured energy losses of the charged particles traversing the 
77-ring divided by the average energy loss per charged par- 
ticle. The charged particle pseudorapidity distribution is 
corrected for (5-electrons produced by lead ions traversing 
the 213 /im thick 208p|^ target foil. On average these elec- 
trons give rise to roughly 11 additional hits in the SPMD. 
This contribution has been determined from beam trig- 
gers where no inelastic interaction took place. The sys- 
tematic error of dNch/drj relates to the uncertainty in the 
determination of the total energy loss of the charged par- 
ticles in the SPMD and to the correction for (5-electrons. 
The uncertainty in the energy loss measurement is esti- 
mated to result in a 3% centrality independent system- 
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atic error. The correction for (5-electrons at midrapidity 
{dNs/drj\mid ~ 9) is assumed to be known with an ac- 
curacy of 10% and contributes significantly to the total 
uncertainty only for peripheral reactions. However, the 
dNch/dr] distributions as shown in figure |6| exhibit a slight 
asymmetry around midrapidity. We can force the dNch/dr/ 
distribution to be symmetric by arbitrarily increasing the 
subtracted (5-electron contribution. For central events the 
^-electron yield has to be increased by roughly 80%. This 
factor decreases when going to semi-central and periph- 
eral events. It's difficult to imagine a physical reason for 
a much larger (5-electron production than what was mea- 
sured in target-out events. However, by making these ex- 
treme assumptions we estimate the centrality dependent 
error of dNch/dr] by comparing the dNch/dr] distributions 
with (5-electrons subtracted as measured in beam events 
with the dNch/dr] distributions that were forced to be 
symmetric. For peripheral and semi-central events this un- 
certainty typically is of the order of 3—4% and decreases to 
2% in central events. The correction of dNch/dr] measured 
with the SPMD in peripheral events due to interactions 
outside the target is typically of the order of 15%. The un- 
certainty of dNch/dr] due to this correction is estimated 
to be around 3%. 

The photon distributions used in this analysis were 
measured with the LEDA spectrometer in the rapidity 
interval 2.3 < 77 < 3.0. This detector is located 21.5 m 
from the target and consists of 10080 leadglass modules 
each read out by a photomultiplier. A streamer tube array 
placed directly in front of LEDA was used as a charged 
particle veto detector to correct for the charged hadron 
contamination in the leadglass. The remaining correction 
for neutrons and anti-neutrons has been made based on 



simulation results using the GEANT package 1 15 . The de- 
tection efficiency of photons in LEDA is based on GEANT 
simulations and experimental data in order to take into ac- 
count the effects of overlapping showers which can result 
in a shift in the measured transverse momentum. These 
corrections require high statistics and therefore only 8 cen- 
trality classes have been used for the photon analysis here. 
These are the same 8 centrality classes used in the analy- 
sis of the scaling of neutral pion production presented in 
1^. The systematic error on the photon and neutral pion 
multiplicities is estimated to be « 10% mainly originating 
from corrections for efficiency and contamination. 



3 Model Calculations 

In the present analysis, the photon and charged particle 
scaling has been investigated with the centrality of the 
Pb-|-Pb collision determined from the transverse energy, 
measured with the MIRAC calorimeter. However, in 
the Et scaling analysis, the forward energy Ep of projec- 
tile spectators measured with the Zero Degree Calorimeter 
has been used for the centrality selection, in order to avoid 
auto-correlations. Twenty-one centrality classes have been 
defined based on the measured Et- Each class corresponds 
to 5% of the minimum bias cross section, with an addi- 
tional very central class corresponding to the 1% most cen- 



tral events. The ZDC cannot resolve the centrality of very 
peripheral collisions as well as the MIRAC calorimeter and 
therefore the Et scaling analysis is limited to centrality 
classes which correspond to more than approximately 40 
participants. 

The number of participants Npart and collisions Ncoii 
for a given centrality class have been determined from 
a simulation based on the event generator VENUS 4.12 
Control calculations were made with the FRITIOF 
event generator [ p9| . The definition of Npart and Ncou in 
these models is based on a geometrical picture (Glauber 
model) where nucleons travel on straight-line trajectories 
and the nucleon-nucleon cross section is independent of 
the number of collisions a nucleon has undergone before. 
The approach of using event generators allowed to take 
into account the energy resolution of the ZDC and MIRAC 
calorimeters. Furthermore, the minimum bias trigger effi- 
ciency has been included in the simulation. Figure y shows 
a comparison of the measured Et and Ep distributions 
to the VENUS simulations. The overall agreement in the 
Et distribution between the data and the model is good. 
However, the VENUS prediction extends to slightly higher 
transverse energy for the most central reactions. In the 
forward energy distribution, the strong peak for periph- 
eral reactions is not precisely reproduced, while the gen- 
eral shape is quite similar. Also the event-by-event anti- 
correlation of Et and Ep observed in the experimental 
data is in good agreement with the model calculations 
(see figure ^. 

To obtain a robust estimate of Npart and Ncoih which 
is less sensitive to discrepancies in the energy distribu- 
tions, the centrality classes in the model have been chosen 
to represent the same absolute cross section as the data. 
The effect of the centrality cuts on the distributions of the 
number of participants can be seen from figure ^, where 
distributions of Npart for centrality classes corresponding 
to fractions of the minimum bias cross section of 0-1%, 1- 
5%, 5-10%, 10-20%, 20-40%, 40-60%, 60-80% and 80-100% 
are shown. One can see that the limited acceptance, the 
detector resolution, and the fluctuations in particle pro- 
duction as implemented in the model lead to an overlap of 
the distributions of adjacent centrality classes. Neverthe- 
less, it is observed that even a strong cut on the 1% most 
central reactions yields a significantly different selection 
than e.g. the 5% most central reactions. 

Figure || shows a summary of the RMS-values versus 
the average values of the distributions of Npart for the 
centrality selections as used in the later analysis in this 
paper. While the classes selected with Ep have a slightly 
smaller width for very central reactions, for peripheral re- 
actions the resolution of the selection is much better using 
Et- 

The precise number of participants or number of col- 
lisions may, however, depend on the specific model as- 
sumptions. We have therefore performed a detailed study 
of the influence of these assumptions and other sources of 
systematic errors for the analysis presented here. A sum- 
mary of these studies is given in appendix 
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Fig. 2. Distributions of tlie transverse energy Et as measured 
in MIRAC (upper graph) and the forward energy Ef as mea- 
sured with the ZDC (lower graph) in 158-A GeV Pb+Pb col- 
hsions. Predictions of the event generator VENUS 4.12 are 
included. 



4 Results 

The pseudorapidity distributions for the transverse en- 
ergy, the charged particle multiplicity, and the photon 
multiplicity are shown in figure ^ for five centrality classes. 
All distributions have been corrected for possible contri- 
butions of reactions upstream and downstream of the tar- 
get and in the target frame by subtracting the respective 
yield determined in target-out runs, see for further 
details. To obtain dEx/drjlmax the measured data points 
have been reflected at midrapidity (77cm = 2.91) and fitted 
with a Gaussian. 

A first impression of the centrality dependence of Et 
and charged particle production can be obtained by nor- 
malizing the yields to the number of participants. This is 
shown in figure 0. Both for Et and the charged multiplic- 
ity the yield per participant increases when going from 
peripheral to more central reactions. As described in the 
appendix we have also used the FRITIOF model instead 
of VENUS to calculate the number of participants and 
the number of nucleon-nucleon collisions. The results of 
these two calculations for the number of participants are 
compared in the appendix (fig. nw and used to obtain the 
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Fig. 3. Event-by-event distributions of the transverse energy 
Et vs. the forward energy Ep in 158-^4 GeV Pb-|-Pb collisions. 
The left plot shows the experimental data and the right plot 
predictions of the event generator VENUS 4.12. 
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Fig. 4. Distributions of the number of participants Npart ob- 
tained from VENUS 4.12 for different centrality classes selected 
by the transverse energy Et- 



estimated uncertainty of the number of participants which 
was used in the calculation of the error bars in figure ^ 

The scaling behavior is studied in more detail in fig- 
ure js] which shows the dependence of the Et and Nch 
pseudorapidity densities at midrapidity on the number of 
participants. The scaling behavior of these observables was 
parameterized as 



dX 
drj 



NparuK^U, X^ET,N,h. 



mid 



(1) 



This functional dependence gives a reasonable descrip- 
tion of the data for the entire centrality range. Taking 
the number of participants from the VENUS calculation 
(denoted as calculation A in the appendix) the charged 
particle scaling can be described by a scaling exponent 
Up = 1.08. The calculations of the number of participants 
using VENUS, or VENUS with an experimentally deter- 
mined nucleon density distribution, or FRITIOF (denoted 
as calculations A, B, and F in the appendix) arc all based 
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Fig. 5. RMS-values of the number of participants Npart ob- 
tained from VENUS 4.12 for different centrality classes selected 
by the transverse energy Et and the forward energy Ef- 
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Fig. 6. Pseudorapidity distributions of transverse energy, 
charged particles and photons measured in 158-^4 GeV Pb-|-Pb 
collisions of different centrality. Photons were measured above 
a lower transverse momentum threshold of 180 MeV/c. The 
open symbols have been obtained by reflecting the measured 
data points at midrapidity. 



on reasonable assumptions. These three calculations give 
slightly different participant numbers. We quote the aver- 
age of the corresponding three scaling exponents as the 
final result. As described in the appendix the estimated 
uncertainty of due to the uncertainty of the number of 
participants is 0.036. In the fit from which we obtain the 
value of Qfp we only take the statistical error of the data 
points into account. Since the statistical errors are small 
the fit error is negligible. In order to estimate the influ- 
ence of the centrality dependent errors of dNch/drj on ap 




100 200 300 400 

N , 



Fig. 7. Transverse energy and charged particle yields in 
158-^ GeV Pb-|-Pb reactions normalized to the number of par- 
ticipants as a function of the number of participants. For the 
number of participants an uncertainty as indicated in figure 
was assumed. Furthermore, all centrality dependent errors of 
dEr/drj and dNch/dr] as described in section ^ were taken 
into account in the calculation of the error bars. The brackets 
indicate the centrality independent uncertainty of the overall 
dEx/drj and dNch/drj scale. 



we systematically move the data points within the error 
bars and repeat the fit. From this procedure we estimate 
an uncertainty of 0.02 for ap due to the centrality depen- 
dent errors of dNch/drj. Adding all errors in quadrature 
we finally obtain ap = 1.07 ±0.04 for the charged particle 
scaling. 

For the Et scaling we perform a similar error analysis. 
An additional uncertainty comes from the assumed cen- 
troid of the cLEt / drj distribution. Primarily due to massive 
particles like protons and neutrons, the difference between 
pseudorapidity and rapidity could lead to an increase of 
the centroid relative to midrapidity [ijmid = 2.91). By 
varying the assumed dEx/drj centroid position in the r]- 
range 2.91 ± 0.3 a corresponding error of 0.02 was esti- 
mated for the Et scaling exponent ap. Adding all errors 
in quadrature we obtain ap = 1.08 ± 0.06. 

In more detail, the relative scaling for different cen- 
tralities can be judged from the local scaling exponent 
diocai shown on the right hand side of figure ^. These 
have been obtained from a fit of 5 neighboring Et and 
Nch data points. For the charged particles it can be seen 
that the scaling remains approximately constant over the 
whole centrality range. For Et the local scaling exponent 
appears to be almost constant in the range Npart > 100. 
However, below Npart ~ 100 the local scaling exponent 
seems to increase slightly. 

Considering the scaling with the number of binary 
nucleon-nucleon collisions and averaging the results deter- 
mined with calculations A, B, and F we obtain a similarly 
good description with ac = 0.82±0.03 for the charged par- 
ticle scaling. For the Et scaling we obtain ac — 0.83±0.05. 
These results are not surprising since for symmetric sys- 
tems one naively expects a scaling relation to hold between 
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Fig. 8. Pseudorapidity density of Et and Nch at midrapid- 
ity as a function of the number of participants. The partici- 
pant numbers shown here were calculated in a VENUS simu- 
lation (denoted as calculation A in the appendix). Both for 
Et and N^h the fit results obtained from a fit of all data 
points are shown. In order to demonstrate the stronger than 
linear increase of the data points a linear extrapolation of the 
charged particle multiplicity in peripheral Pb-|-Pb reactions 
{Npart ~ 50) is shown as a dotted line in the lower left plot. 
The scaling behavior can be verified in more detail on the right 
panel where the local scaling exponents are shown. The local 
scaling exponents have been obtained from a fit of 5 neighbor- 
ing Et and Nch data points. 

the number of collisions and the number of partieipants: 
Ncoii oc Np^^^ with Qfcp = ^. (2) 

Fits of the parameters extracted from VENUS simulations 
indeed yield a value of Ucp = 1.28 which is close to the 
above value. For a scaling with A^p^^ this would lead to a 
behavior as N^^^j^- 

A similar analysis can be performed on the data ob- 
tained from the VENUS simulation itself. The results of 
such an analysis are displayed in figure ^. It is observed 
that the scaling exponents are higher in the simulation, 
and that both for Et and Nch the exponent a shows a 
tendency to increase with centrality. We note here that the 
extraction of the number of participants and the number 
of binary nucleon-nucleon collisions is almost completely 
independent of the scaling exponent present in the under- 
lying event generator. As described in section ^ this is due 
to the fact that the centrality classes defined in the model 
calculations correspond to the same absolute cross section 
as the experimental centrality classes. 

It is interesting to extrapolate this scaling towards 
smaller system sizes and compare to the expectation from 
pp collisions |^ . This has been done in figure It can 
be seen that the scaling obtained from charged particle 
production in Pb-|-Pb collisions extrapolates nicely to pp 
collisions. In particular, there is no threshold effect visible 




400 



Fig. 9. Pseudorapidity density of Et and Nch at midrapidity 
as a function of the number of participants as in figure pi from 
VENUS 4.12 simulations. 

in the charged particle multiplicity when going to central 
Pb+Pb collisions. 




Fig. 10. Pseudorapidity density of Nch at midrapidity as a 
function of the number of participants for p-l-p and Pb-|-Pb 
collisions. The number of participants shown here were calcu- 
lated using VENUS (calculation A in the appendix). The fit 
function plotted as a solid line is the same as obtained in fig- 
ure ^. The dashed function is the fit result using the participant 
values from the FRITIOF calculation (calculation F). On the 
right hand side the ratio of the data (using participants from 
VENUS) to the fit function is shown. 

In a recent publication we have discussed the sys- 
tematics of inclusive photon production as measured with 
the WA98 Photon Multiphcity Detector. There it was 
found that the pseudorapidity density of photons at midra- 
pidity scales with the number of participants as A'^J^^j^^'^^. 
This is slightly larger, but consistent with the exponent 
from the present analysis of Et and Net- Photon produc- 
tion can be investigated in further detail with the photon 
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Fig. 11. Exponents for a scaling with the number of partici- 
pants and nucleon-nucleon colhsions for photons measured in 
LEDA as a function of the transverse momentum. Only reac- 
tions with Npart > 30 have been used. The data points shown 
were obtained using the number of participants from VENUS. 
The error bars reflect the uncertainty of the photon measure- 
ment, the uncertainty in the calculation of Npart/ N^oii and the 
fit error. 



spectrometer LEDA. It can very naturally be studied as a 
function of the transverse momentum. The corresponding 
scaling exponents a extracted are shovi^n in figure nil The 
photon measurement in LEDA suffers from larger system- 
atic uncertainties at low momenta, so only photons for 
Pt > 500 MeV/c have been considered. At a transverse 
momentum oipT ~ 500 MeV/c the inclusive photon yield 
shows a scaling behavior similar to that observed for Et 
and Nch- However, the extracted scaling exponents tend 
to rise with increasing and at pt ~ 2 GeV/c the scaling 
can be described as ^ N^art ^'^d ~ N^m. 

Since a large fraction of the inclusive photons origi- 
nates from the decay of neutral pions it is of interest to 
compare the scaling of photons to that of neutral pions 
which have already been discussed in [|| . These data have 
been reanalyzed and the results of the scaling expo- 
nents with respect to the number of participants are shown 
in figure The values are nearly constant at a k, 1.1 
with a tendency to decrease towards higher transverse mo- 
menta. It should be noted that the extracted exponents 
are smaller compared to the values given in This is 
mostly due to a more sophisticated calculation of the num- 
ber of participants used here. 

We have recently published results on the production 
of direct photons in 158- A GeV Pb-fPb colhsions ||2^,|l|. 
A significant yield of direct photons at pt > 1-5 GeV/c 
is observed in central collisions, while in peripheral colli- 
sions the photon production is consistent with the yield 
expected from the decays of neutral pions, -q mesons, and 
other hadrons. At pr ~ 2 GeV/c the direct photon yield 
in central collisions amounts to roughly 20% of the pho- 
tons from hadronic decays. The scaling of the neutral pion 
yield presented in figure |l^ appears to be consistent with 
the scaling of charged particles and the transverse en- 
ergy, while the photon yield at higher transverse momenta 



seems to increase more strongly with the number of par- 
ticipants. We assume that for a relatively small yield of 
direct photons in central collisions, the central! ty depen- 
dence of the inclusive photons can still be described with 
a scaling law as in equation It's then obvious that the 
production of direct photons in central Pb-|-Pb collisions 
necessarily increases the scaling exponent a of the inclu- 
sive photons at high transverse momenta. In this respect 
the behavior of the extracted scaling exponents for inclu- 
sive photon and neutral pion production are consistent 
with the direct photon excess observed in central Pb-t-Pb 
collisions. However, due to the uncertainties of the scaling 
exponents for photons and neutral pions it is not possi- 
ble to draw quantitative conclusions about direct photon 
production from figures |ll] and [ij. 
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Fig. 12. Exponents for a scaling with the number of partic- 
ipants and nucleon-nucleon collisions for neutral pions mea- 
sured in LEDA as a function of the transverse momentum. 
Only reactions with Npart > 30 have been used. The data 
points shown were obtained using the number of participants 
from VENUS. The error bars refiect the uncertainty of the tt" 
measurement, the uncertainty in the calculation of Npart/ Ncoii 
and the fit error. 



In simple multiple collision models a heavy ion reaction 
is regarded as a sequence of independent nucleon-nucleon 
collisions which can be described as in free space |^ , ^ . 
After a projectile nucleon suffers an inelastic collision the 
assumption of local baryon number conservation assures 
that a baryon-like object is still present. This baryon-like 
object is assumed to contribute to the particle produc- 
tion in subsequent collisions with the same cross section 
as the initial nucleon. In this picture the contribution of 
each nucleon-nucleon collision is added incoherently which 
leads to a linear scaling of Et and particle production with 
the number of binary collisions. If the energy degradation 
in each nucleon-nucleon reaction is taken into account a 
reasonable description of Et and particle production can 
be obtained. The approximate scaling as N^^f for the 
transverse energy and charged particles may be used to 
obtain information on the average energy degradation in 
a nucleon-nucleon collision. 
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Fig. 13. Pseudorapidity density of Nch at midrapidity in 
Pb+Pb collisions normalized to the number of collisions and 
the charged particle density in pp collisions. The solid line 
shows a fit with equation (Ml). 

A simple way of investigating this hypothesis is to 
study the particle production per binary collision as a 
function of the effective thickness x of the two nuclei, 
which might be characterized in the following way |23,E4,E5[ 



-{AA) = Ncou ■ f[x) ■ — — (pp). 



dr] 



dr] 



(3) 



We have attempted to fit the data in figure 
the form: 



Here f{x) should describe the effect of energy degrada- 
tion on particle production with x being a suitable thick- 
ness variable ~ we have chosen x = Ncoii/Npart- Figure ^ 
shows the pseudorapidity density of Nch at midrapidity in 
Pb-|-Pb collisions normalized to the number of collisions 
and the charged particle density in pp collisions. The data 
show a continuous decrease of the multiplicity per colli- 
sion for increasing x, i.e. the more collisions a participant 
suffers, the smaller is the contribution of each collision 
to particle production. As an example, this can be illus- 
trated more quantitatively for the centrality classes with 
Ncou/Npart ~ 1, i-e. for the case that each participating 
nucleon on the average suffers two nucleon-nucleon colli- 
sions. For these reactions the actual charged particle den- 
sity at midrapidity is 40% lower than one would expect 
from a linear scaling of the charged particle density in pp 
reactions with the number of binary nucleon-nucleon col- 
lisions. In a simple multiple collision picture this means 
that on the average the second collisions of each partici- 
pant contributes only 20% of the yield of the first collision 
to the charged multiplicity at midrapidity. 

In proton-nucleus reactions multiplication of the charged 
particle multiplicity observed in pp reactions by the num- 
ber of binary nucleon-nucleon collisions overestimates the 
measured multiplicity by 20—30% p^ , |27t . In case of proton- 
nucleus reactions the target participants suffer exactly one 
nucleon-nucleon collision. This is of course not true in AA 
collisions. Nevertheless, it is interesting to apply the same 
recipe to heavy ion reactions. Figure 03 shows that mul- 
tiplying the pp yield with the number of binary nucleon- 
nucleon collisions gives as much as 60% too many charged 
particles in central Pb-|-Pb reactions. 
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(4) 



This function basically represents a first order Taylor ap- 
proximation of the function 1// in equation (^. A good 
fit is obtained with a = 1.19 ± 0.05. 

In the case of the Wounded Nucleon Model the even- 
tual fragmentation of an excited nucleon after an inelastic 
collision is not affected by further collisions with other nu- 
cleons, no matter how many times it is successively struck. 
The observation of a stronger than linear increase of Et 
and Nch with the number of participants indicates that 
this model is only approximately correct. It can be clearly 
seen, regarding the two possible centrality variables in- 
vestigated here, that the number of participants is better 
suited, because the scaling exponent is closer to one com- 
pared to the number of collisions. 

With the present data it is possible to determine the 
average transverse energy per charged particle at midra- 
pidity 

{ET)/{Nch)\rmd = (dET/dr]\,md)/{dNch/d7]\,md), (5) 

a quantity that can be seen as a measure of the global 
mean transverse momentum averaged over all particle species. 
{Et}/ {Nch)\mid is plotted in figure |lj as a function of cen- 
trality, represented by the number of participants. For this 
figure both dEx / dr] and dNch/drj were evaluated in iden- 
tical centrality classes, defined with the forward energy 
Ep. 
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Fig. 14. {Et) / {Nch)\mid as a function of the number of par- 
ticipants. For comparison results of VENUS 4.12 calculations 
are included. The error bars indicate the centrality depen- 
dent errors of the experimental result. The overall uncertainty 
of {Et) / {N ch)\mid which is mainly due to the uncertainty in 
{dET / drfjmid is indicated by brackets. 



{Et) I {Nch)\mid appears to increase up to a system size 
of Npart ~ 100 which corresponds to an impact parameter 
of 6 w 9 fm. For more central collisions {Et) / {N ch)\m,id 
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levels off at a value of 0.80 GeV. This value is slightly 
higher than the maximum {Et) / {N ch)\mid ~ 0.66 GeV 
observed in 200-^ GeV S+Au and S+Al reactions 0. 
VENUS 4.12 predicts a qualitatively similar behavior for 
{Et) I (Nch) \mid, while the absolute value is approximately 
100 MeV lower at Npart ~ 100. One may also note that 
the VENUS results continue to rise by « 50 MeV when 
going from Npart ~ 100 to Npart ~ 400, while the exper- 
imental data appear to be completely flat in this region. 
A similar saturation with increasing number of partici- 
pants as observed here has been seen in the (truncated) 
mean px of neutral pions with px > 400 MeV/c produced 
in Pb-I-Pb reactions |^. A natural explanation of such a 
behavior would be the assumption that thermalization is 
reached once the system exceeds a certain minimum size. 



5 Conclusions 

We have analyzed the dependence of transverse en- 
ergy and charged particle pseudorapidity distributions 
and photon transverse momentum spectra in 158- A GeV 
Pb-|-Pb collisions on the number of participants and 
the number of binary nucleon-nucleon collisions. A scal- 
ing behavior as N^aZ^°'°'^ and A^°„ff describes the 
charged particle production over the whole impact param- 
eter range. The Et production was studied for collisions 
with more than approximately 40 participants. In this cen- 
trality range the Et production scales as N^^^^^^'^^ and 



0.83±0.05 

coll 



n: 

Photons at pT ~ 500 MeV/c show a scaling behavior 
similar to the scaling of Et and Nch- The pr-dependence 
of the photon scaling was studied and a rise of the ex- 
tracted scaling exponents with increasing transverse mo- 
mentum was found. 

We have studied the transverse energy per charged 
particle as a function of centrality and found an indication 
of an increase from peripheral to semi-central collisions 
with approximately 100 participants with a subsequent 
saturation for larger systems. 

While the global variables like Et and charged particle 
multiplicity seem to scale closer to the number of partic- 
ipants than to the number of binary collisions, there is a 
clear participant scaling violation compared to a purely 
linear dependence. This scaling violation might e.g. have 
consequences for the suppression of J/ip production from 
comovers, since the implied central particle densities are 
considerably larger than estimated based on a linear scal- 
ing. 



A Systematic uncertainties in model 
calculations 

In order to obtain an estimate of the systematic uncer- 
tainties in the calculation of the number of participants 
and the number of nucleon-nucleon collisions we have var- 
ied several assumptions in the model calculations. The 
dependence of the particle and transverse energy yield on 



the number of participants and collisions is described with 
the scaling exponents ap and ac in this paper. In this sec- 
tion we investigate how the extracted scaling exponents 
ttp and ac for the charged particle yield (with similar con- 
clusions for the other observables) are affected by the dif- 
ferent model assumptions. 

In the VENUS 4.12 simulations used to obtain the 
number of participants and collisions we have varied 

— the parameterization of the nucleon density distribu- 
tion, 

— the energy resolution of MIRAC and 

— the minimum bias cross section. 

As a cross check we have also calculated the number of par- 
ticipants and collisions using the event generator FRITIOF 

The nuclear density profile used in VENUS and FRITIOF 
is an effective parameterization using a Woods-Saxon shape: 



p{r) 



I 



Po 



H-exp(i^) 



(6) 



However, the two models make slightly different assump- 
tions for the nuclear radius R and the diffuseness param- 
eter a. The VENUS parameterization is 



i?VEN = (l.lMi/^ 
avEN = 0.54 fm 



1.61^-1/3) 



(7) 



which results in i?vEN ~ 6.78 fm for a lead nucleus. The 
radius parameter in FRITIOF for nuclei with ^ > 16 is 
calculated as 



i?pj^j = ro ■ A^/^ with 

ro = 1.16 • (1 - 1.16^-2/3^) fj^^ 



(8) 



This gives i?FRi « 6.65 fm for lead nuclei. The diffuse- 
ness parameter a is taken to be slightly yl-dependent in 
FRITIOF and hcs in the range 0.47 fm - 0.55 fm. For lead 
nuclei FRITIOF uses 



flpRi — 0.545 fm. 



(9) 



Electron scattering experiments have shown, however, 
that the density distribution has a slightly more compli- 
cated structure Q For a comparison we will use a parame- 
terization fitted to electron scattering data on ^"^Pb pre- 
sented in lESi: 



Poxp(?') = Pa 



Ci -I- C2r + c^r'' 



1 -I- exp 



\ flex 



(10) 



with 



ci = 0.0633, C2 = -0.002045 fm-\c3 = 0.000566 fm'^ 
and 



6.413 fm, flc 



0.5831 fm. 



^ Strictly this applies only to the charge distribution. The 
true nucleon distribution, especially in the inner part of the 
nuclei, is not experimentally accessible. 



11 



The three different distributions are shown in figure 15 , It 
can be seen that the overall agreement is quite good. The 
default distribution of VENUS has a slightly larger radius, 
while the parameterization of the experimental data shows 
small oscillations compared to the other two distributions. 




Fig. 15. Different nuclear density distributions used in the 
calculations of the number of participants and the number of 
collisions. The solid line shows the distribution implemented in 
the VENUS 4.12 Monte-Carlo model (equation (||) with the pa- 
rameters (0) ) , the dashed line shows a parameterization of the 
charge distribution obtained from electron scattering (equa- 
tion ([lo|)) and the dotted line shows the density used in the 
FRITIOF model (equation with the parameters (^ and 



With respect to the experimental resolution we have 
varied the energy resolution of the MIRAC calorimeter 
in the simulations. The measured values of the resolution 
Ip^ are for the electromagnetic section: 



tlcm ^ 17.9% 

and for the hadronic section: 



fhad 



46.1% 



This has been arbitrarily worsened to 



.,had 



85% 



E y/E/GeV 



(11) 



(12) 



(13) 



for both sections of the calorimeter. 

For an accurate determination of the number of partic- 
ipants and the number of nucleon-nucleon collisions it is 
necessary that the experimental minimum bias threshold 
is reproduced in the simulation. In section |^ we stated the 
two sources in the determination of the experimental min- 
imum bias cross-section: the error of the target thickness 
and the error due to subtraction of interactions outside the 



target. In our approach of calculating the number of par- 
ticipants and nucleon-nucleon collisions we define central- 
ity classes based on the simulated Et (in case of the Nch 
scaling analysis) which correspond to the same absolute 
cross-sections as the respective classes for the measured 
Et- The uncertainty of the minimum bias cross-section 
due to the error of the target thickness therefore directly 
leads to an error of the iVpart-values which affects the 
entire centrality range. However, the uncertainty in the 
correction for interactions outside the target affects the 
experimental da /cLEt distribution only in the peripheral 
range. More precisely, since the maximum Et measured in 
target-out events is around 75 GeV, according to table |l| 
only the range Npart < 50 is affected. 

The measured target thickness is 213 ± 3 fim. The un- 
certainty of the target thickness contributes a relative er- 
ror of roughly 1.5% to the error of the experimental min- 
imum bias cross section. In order to check the influence 
on the extracted scaling exponent we arbitrarily increase 
the minimum bias threshold in the simulation by roughly 
2% to CT^™ = 6386 mb. We have furthermore checked the 
error of ap due to the uncertainty of the minimum bias 
cross section that relates to the subtraction of non-target 
contributions. To do this, we have used an Et distribu- 
tion without correction for reactions outside the target 
to translate the experimental Et cuts into cross section 
cuts. The apparent minimum bias cross section for this 
Et distribution was cr„ifc = 6530 mb. 

The different calculations of the number of participants 
and collisions are summarized in the following list: 

A VENUS 4.12 calculations using the standard settings 
for the density distribution, the experimental resolu- 
tion and the minimum bias cross section. 

B VENUS 4.12 calculations as in case ^ with a modified 
density distribution according to equation ( p^ ) . 

C VENUS 4.12 calculations as in case ^ with a modified 
MIRAC resolution. 

D VENUS 4.12 calculations as in case|^ with an a mini- 
mum bias cross section increased to a^™ = 6386 mb. 

E VENUS 4.12 calculations with cross section cuts de- 
rived from an experimental i^T-spectrum that was not 
corrected for interactions outside the target (apparent 
minimum bias cross section: amb = 6530 mb). 

F FRITIOF calculations using standard settings as in 
calculation A. 

As an example, the number of participants and col- 
lisions from calculation A and F are compared in figure 
|l6| . In peripheral Pb-fPb reactions VENUS gives up to 
10% more participants than FRITIOF whereas in central 
reactions both simulations yield almost identical results. 
Almost independent of centrality the number of collisions 
from VENUS is roughly 10% higher than the FRITIOF 
result. The results of calculation A and F together with 
the experimental Et intervals are given in table ^ 

For the case of the Nch scaling the impact of the differ- 
ent model assumptions on the extracted scaling exponents 
a is summarized in table ^. Considering the scaling with 
the number of participants we take the average of the ap 
values from calculations A, B and F as our final result 
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Fig. 16. Relative difference of tfie number of participants (left) 
and the number of nucleon-nucleon collisions (right) calculated 
with VENUS (calculation A) and FRITIOF (calculation F). 
The dashed line on the left plot indicates the assumed uncer- 
tainty of the participant scale that is used in the calculation of 
the error bars in figure ^ 



% 

of c.s. 


(GeV) 


^ part 1 coll 

VENUS 


-^part 1 coll 

FRITIOF 


1 


398.8 


380.7 


810.7 


380.3 


717.4 


5 


355.8 


355.8 


739.4 


355.4 


662.3 


10 


313.1 


310.9 


621.7 


313.0 


562.9 


15 


275.2 


269.7 


518.5 


270.0 


469.9 


20 


239.8 


233.3 


429.5 


232.8 


389.4 


25 


208.0 


202.0 


357.2 


201.0 


324.5 


30 


179.2 


174.2 


293.6 


172.3 


264.7 


35 


153.7 


149.3 


240.1 


147.3 


218.6 


40 


130.3 


128.2 


197.3 


125.7 


177.2 


45 


109.7 


109.3 


159.4 


105.4 


141.8 


50 


91.2 


91.4 


126.1 


89.2 


114.6 


55 


74.8 


76.2 


99.2 


73.3 


89.2 


60 


60.4 


62.6 


76.7 


60.3 


69.1 


65 


47.9 


51.2 


59.1 


49.2 


53.4 


70 


37.0 


41.3 


44.8 


39.1 


39.6 


75 


27.9 


32.4 


32.6 


30.9 


29.8 


80 


20.5 


25.5 


24.2 


23.9 


21.5 


85 


14.7 


19.5 


17.3 


17.8 


15.1 


90 


10.3 


14.6 


12.1 


13.1 


10.4 


95 


6.9 


10.9 


8.6 


9.8 


7.4 


100 


0.0 


8.3 


6.2 


7.5 


5.4 



Table 1. The number of participants and binary collisions for 
different centrality classes obtained with the measured trans- 
verse energy in Pb+Pb collisions calculated from VENUS and 
FRITIOF (see text). 



since all three calculations are based on reasonable as- 
sumptions. The maximum difference of 0.03 between the 
ap values from these calculations is taken as one contribu- 
tion to the systematic error. By adding the deviations of 
the results from calculations C, D and E from the mean 
value 1.07 in quadrature we estimate the total systematic 
error of ap related to the uncertainty of the number of 
participants to be 0.036. The same prescription yields an 



calculation 


Cp 








A 


0.83 


1.08 


1.94 


0.83 


B 


0.88 


1.07 


2.15 


0.81 


C 


0.86 


1.07 


1.99 


0.83 


D 


0.98 


1.05 


2.13 


0.83 


E 


0.87 


1.07 


2.01 


0.82 


F 


0.97 


1.05 


2.13 


0.83 



Table 2. Influence of different model assumptions on the ex- 
tracted exponents Up and Qc which describe the scaling of the 
charged particle yield with Npart and Ncoii according to equa- 
tion hi. In addition to the scaling exponents we quote the pro- 
portionality constants for the scaling with Npart (cp) and Ncou 
(cc) for each calculation. The full centrality range was used in 
the flt of equation to the measured charged particle yields. 



uncertainty of 0.024 for the exponent ac that described 
the scaling with the number of binary nucleon-nucleon 
collisions. 
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